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1 Introduction 

Where does the force of gravity come from? Anyone unfortunate enough to be near 
a strongly gravitating object, such as a black hole, would say "from the curvature 
of spacetime". In a weaker gravitational field, we measure fluctuations about some 
background spacetime p. These fluctuations correspond to spin-2 particles called 
gravitons. Traditionally, we now say that "gravity comes from the exchange of mass- 
less gravitons" . We believe this because it reproduces Newton's Law of gravitation, 
and that is well tested experimentally. But how well is Newton's Law really tested? 
The truth is that our experimental knowledge of gravity only covers distances be- 
tween 0.2mm and lO^^cm. The lower bound is as a result of Cavendish experiments, 
and the upper bound corresponds to 1% of the current Hubble length. In units of 
h = c = 1, this is equivalent to an energy scale 



Perhaps, therefore, we should alter our previous statement. Gravity at the experi- 
mental scale (H} could be due to the exchange of massless and/or massive particles, so 
long as the masses are less than 10~^^ eV. Indeed, massive gravitons appear automat- 
ically in some higher derivative gravity theories 0. The simplest non-trivial scenario 
would contain a single massive graviton PP. If we have a combination of gravitons 
of many different masses, we have multigravity [3 IH El E] • In this paper, we 
will encounter bigravity [HI GDI- This is the simplest example of multigravity, in which 
gravity is mediated by both a massless graviton and a single ultralight graviton. 

Multigravity is clearly of interest from a purely theoretical point of view. However, 
it is also of interest to phenomenologists because it predicts new gravitational physics 
at very large distances. At distances beyond the Compton wavelength of the ultralight 
mode, the ultrahght mode is turned off, and gravity is mediated by the massless mode 
alone. Large distance modifications of gravity have been in vogue recently, as they 
could offer an explanation to the current acceleration of the universe fTJ ^1 El ITU 



Unfortunately, there are number of problems with many existing models of mod- 
ified gravity. One very serious problem is the presence of ghosts (see, for example, 
lainillHlEll^nj). Although the "AdS brane" model in ^ is ghost-free, modifications 
of gravity are hidden behind the AdS horizon. The six dimensional model of Kogan 
et al is also thought to be ghost-free, although this has not been confirmed as 
the model is too difficult to work with. The DGP model meanwhile, is simple, 
ghost-free, and even predicts new infra-red physics, that could one day be observ- 
able. However, it is not a bigravity model; gravity is due to a resonance of continuum 
modes. 

Let us describe the DGP model in more detail. It is given by the following action 



In other words, we have a four- dimensional Minkowski brane embedded in a five- 
dimensional Minkowski bulk. The key ingredient is the brane localised curvature, 



10"^^ eV < p < 10"^ eV 



(1) 



cniiinj- 
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i?(7). This could be generated by quantum corrections, if matter were present on 
the brane. Locahsed curvature can also appear in string theory j22ll2Sl- In the DGP 
model, Newton's Law is reproduced up to a distance r = mpi/2M^. Beyond this scale 
the behaviour of gravity is five- dimensional. 

In this paper, we will consider a generalisation of the DGP model. We will add 
a second brane with localised curvature. We will also allow for a bulk cosmological 
constant, and for the branes to have tension. Since the extra dimension is finite, we 
will get a discrete graviton mass spectrum. Our aim is to ask the following question: 
is it possible to obtain an interesting modified theory of gravity without introducing 
ghosts and tachyons? The answer will be "yes". For certain parameter regions, we 
will discover bigravity models that are tachyon and ghost-free, leading to potentially 
observable new physics in the infra-red. It is interesting to note that a single DGP- 
like brane in a compact extra dimension also has a discrete mass spectrum, but does 
not exhibit bigravity |24j . 

The rest of this paper will be organised as follows: in section |2l we will describe 
our set up in more detail, giving the bulk and boundary equations of motion. In 
section |2l we will derive solutions for the background spacetime. We will perturb 
about this background in section EJ arriving at the linearised equations of motion. In 
section El we will focus on the radion mode. This corresponds to fluctuations in the 
brane separation. We will calculate its effective action to quadratic order, in order 
to determine whether or not it is a ghost. For de Sitter branes we will find that the 
radion is tachyonic. In section we will turn our attention to the gravitons. We will 
derive the mass spectrum, and show that we can have observable bigravity. We will 
also check for ghosts by calculating the graviton effective action. In section [3 we will 
play around with our parameters until we find a bigravity model that is ghost-free, 
and might one day lead to observable, new, infra-red physics. Section |H1 will contain 
some final remarks. In particular, we will discuss some important issues, such as the 
famous van Dam-Veltman-Zakharov (vDVZ) discontinuity [53l2n]! and the recently 
discovered strong coupling problems in massive gravity ,27,i;28j. Finally, the appendix 
contains a detailed calculation of the mass spectrum for AdS branes. 



2 The Set Up 

Consider two 3-branes embedded in a five- dimensional bulk spacetime. For simplicity 
we assume that there is Z2 symmetry across each brane. This means that we can 
focus on the manifold, A4, which has one brane as its "left-hand" boundary, dj^i, 
and the other brane as its "right-hand" boundary, dM.R. The entire bulk is now 
given by two copies of M.. 

Our set up is described by the following action 

S = •S'bulk + 'S'brane) (3) 

where the contribution from the bulk is given by 

^buik = 2M3/" ^(i?-2A)+ V [ v^ir« + 2/ (4) 

Jm i=L,R JdMi Jm 
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Here M is the five-dimensional Planck mass, and gab is the bulk metric with corre- 
sponding Ricci scalar, R. We have expilicitly included a bulk cosmological constant, 
A, which can be positive, negative or zero. Any additional contributions from matter 
in the bulk appear in Lm- For i = L, R, g^^ is the induced metric on dAii, and K^ij 
is the extrinsic curvature^ of dAii in A^. Notice that S'buik contains an overall factor 
of two because of the Z2 symmetry. 

The brane part of the action is given by 



^1 



branc 



=L,R "^'9-^' i=L,R 



(5) 



where for i = L, R, ai is the tension of the brane at dAii. Cm gives the contribution 
from any additional matter on the brane. We have also allowed for some localised 
curvature by including the four-dimensional Ricci scalar, TZ^^^ on each brane. This 
is a DGP-like kinetic term |2I], whose contribution is weighted by the distance in 
each case. At this stage we make no assumption about the magnitude, or even the 
sign of the r^. 

The bulk equations of motion are given by the Einstein equations 



R 



■ab 



\Rga 



'^9ab + 



2M3 



- ab 



(6) 



where Tab is the energy momentum tensor for additional matter in the bulk. We also 
need to impose boundary conditions at the branes. For i = L, R, these are given by 
the Z2 symmetric Israel equations [29j at dAii 



K 



(i) 
ab 



a 



12M3 



9ab 



n 



ab a 9ab 



+ 



2M3 



^ab g-' "dab 



(7) 



Here corresponds to the energy momentum tensor for additional matter on dAii, 
and T*^*) is its trace. 



3 Background solutions 

In this section we will derive the metric, cjab, for the background spacetimes. These 
correspond to solutions of the equations of motion when no additional matter is 
present in the bulk, or on the branes. In other words 

Tab = rj) = (8) 

Let us introduce coordinates x° = (x^, z), and assume that the brane at dA4i satisfies 
z = Zi, where zl = and zr = I > 0. In order for us to trust our analysis, we need 
to assume that 

(«' 

1 (i) 1 (i) 

We will adopt the convention that if^j, = ^Cng^^f,, the Lie derivative of the induced metric with 
respect to the outward pointing normal. 
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Now seek solutions of the form 

ds^ = gabdx"'dx^ = {z)g ^^jdx^ dx^ + dz^. (10) 

where g^j^u describes a maximally symmetric four- dimensional spacetime, with Rie- 
mann tensor 

Rfiua/s = {gf_iagui3 — g^i/sgua) ■ (11) 

where k = 0, ±1 for Minkowski, de Sitter and anti-de Sitter space respectively. 

The Einstein equations in the bulk give the following differential equations for 
a{z) 

'a'y_ A2 A a" _ A 
a / 6 ' a 6 



By solving these, we find that for — / < 2; < / [30^ ,211 

^cos [±A;(|z| - c)] A = 6P, k = 1 only 

a{z) = {X[c±k\zW A = 0, /t = 0, 1 (13) 

\. ^g±fc(|^|-c) _ ^gTfc(kl-c)j A = -6A;^ K = 0, ±1, 

where c is an integration constant. We are free to define a length scale by setting 
a(0) = 1. This fixes a relationship between A and c: 

{k/ cos{kc) A = 6A;^, k = 1 only 

1/c A = 0, K = 0, 1 (14) 

2k/ [e^'"' - /te±^^] A = -6P, ^ = 0, ±1, 

For each of the solutions (|T!?jl . the Z2 symmetry about ^ = is explicit, whereas 
the symmetry about z = I can be seen when we identify z with z + 21. When we 
have a negative bulk cosmological constant, these solutions correspond to those found 

in isaEHESHMl. 

The Israel equations on dAii give us the following boundary conditions 

'a{z,) 12M3 'a^{z,) ^ ' 

where 9l = —1 and 9ji = +1. For the special case of A = — 6fc^, k = 0, these 
equations give the usual Randall-Sundrum fine tuning of brane tensions jSZl • Finally, 
we note that the cosmological constant on the brane at dA4i is given by 

A^ 

A, = (16) 

and shall henceforth refer to the branes as flat, de Sitter or anti-de Sitter, when 
K = 0, ±1 respectively. 
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4 Metric Perturbations 



We shall now consider perturbations about the background solutions ()13|) we have 
just derived. We will allow additional matter to be present on the branes, but not in 
the bulk. Let us define gab = Qah + ^Qab to be the perturbed metric. We will work in 
Gaussian Normal coordinates so that 

^9^lz = ^Qzz = 0. (17) 

Since we have no additional bulk matter, we can take the metric to be transverse- 
tracefree in the hulk. In other words, 5g^j_y = Xfiu, where 

v.x): = x'; = o. (18) 

Here V is the covariant derivative with respect to the Minkowski, de Sitter or anti de 
Sitter metric g^y. In this choice of gauge, the linearised bulk equations of motion are 
given by 

^V^ - 4/tA2 2A 



X^u = 0. (19) 



where we have used equation (fT^. 

Unfortunately, we can no longer assume that the branes are fixed at z = Zi. The 
presence of matter on the branes will cause them to bend ^J. In general, they 
will now be positioned at z = Zi + fi{x), for some function /j that depends only on 
the x^. This makes it difficult to apply the Israel equations at the branes. To get 
round this we apply a gauge transformation that fixes dAii, and another that fixes 
dJ^R [Ml, without spoiling the Gaussian Normal condition (fTTj) . This gives rise to 
two coordinate patches that are related by a gauge transformation in the region of 
overlap. We will call the patch that fixes dJ^L, the left-hand coordinate patch, and 
the one that fixes dAi^, the right-hand patch. 

To fix dM.Li 'we make the following coordinate transformation 

z^z-fiix), x^ ^ x^ + a'^{z)d^fL [ -^f-- (20) 



a2(y) 

dJ^L is now fixed at z = 0, although dJ^R \s at z = I + — f^. The metric 
perturbation in this patch is given by 

~ ~ dv a' 

Sg^^u = xit^ = X,u - 2a^V,VjL / ^ + 2-fLg^u (21) 

Jo [y) a 

Similarly, to fix dAi^, let 

z^z-fR{x), x^ ^ x^ + a'^{z)d^fR [ (22) 



a2(y) 

Now we have dM.R at z = I but with OM-l ^t z = fi — Jr. The metric perturbation 
is given by 

dv a' 

59,u = xi^^ = Xm. - 2a'V^VjR / + 2-fR-g,, (23) 

Ji [y) 
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We are now ready to use the Israel equations ((Tj) to give linearised boundary conditions 
at each brane. At dAii, we have 



6i d 



+ 



2 dz 12M3 2a2 



0i + 2r — 
a 



2M3 



h 

(24) 



where everything is evaluated aX z = Zi. In deriving equation (j2H), we have made 
repeated use of equations (fT^ and (fT^ . 

Note that if we take the trace of equation ^I^. and use the fact that Xtiy is 
transverse-tracefree, we get 



0^ + 2r, 



a'{z.j) 
a{zi) 



6M3 



(25) 



This clearly shows that matter on the brane does indeed cause the brane to bend. 

The general perturbation is made up of three parts: a radion, plane wave gravitons, 
and a particular solution that depends explicitly on the matter sources, The 
gravitons will be discussed in section El In the next section we will concentrate on 
the radion. 



5 The Radion 

5.1 Radion wavefunction 

The radion mode, corresponds to fluctuations in the brane separation. It is in 

addition to the fluctuation caused explicitly by the matter on the brane, according to 
equation (j2H)). The radion exists because equation (j^^ can have non-trivial solutions, 
fi, even when T^u = 0. Clearly then, the best way to see the radion is to consider 
a situation where we have no matter on the branes. Let us begin with a transverse- 
tracefree gauge in the bulk 

S9,u = (26) 

with the branes positioned at z = Zi+fi{x). The fields, fi{x), are yet to be determined. 
xt^'^^ satisfies equation (fTIH) . with the following boundary conditions at dA4i 



a' 



(27) 

where, again, everything is evaluated at z = Zi. Motivated by the tensor structure on 
the right-hand side of equation (j^^ . we look for solutions to (fT^ of the form 

x!:;'^ = /^(^) (v^V. + /s:^^^.") (28) 
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where we have introduced a scalar field, which we will call the radion field. Note 

(rad) 



that since xlif is transverse-tracefree, we must have 







(29) 



This means that the radion field, 0, has mass = — 4kA^. When k = 0, the radion 
is massless, whereas for k = —1, it is massive. For k = 1, the radion is tachyonic, 
signalling an instability for de Sitter branes [10] . 

Given the ansatz (j2HI), we solve (fT^ to get an expression for fi{z). 



a + [3z 



aa 



a 



z dy 
(aa'f 



for A = K = 
otherwise 



(30) 



However, these expressions for fi{z) are too general. We find that x^^^'^^ contains 
terms which are pure gauge, and therefore unphysical. In fact, when A = n = 0, 



xt^'^^ can be completely gauged away by the following coordinate tranformation 



z — 



(31) 



When we apply the boundary conditions (j24p . we find that fL = fR = 0. This means 
that there is no physical radion, and no matter-independant brane bending. 

For every other case, the radion cannot be gauged away completely, so we are left 
with a physical radion mode. We can remove the term proportional to a in x^i^i^'^^ by 
taking 

x'' — > x'' + -aaa'df,(j), z ^ z + -anX^cp. (32) 

but we can do nothing about the other term. The physical radion perturbation is 
therefore given by equation ()28|) with 



H{z) 



aa 



dy 



\aa 



(33) 



Notice that we have set /? = 1 by absorbing it into the definition of 0. Finally, the 
boundary conditions (|77|) give the fi in terms of 0(x). Making use of equations (fT^ 
and (jTSJ, we get 



yia{zi)a'{zi) 



dy 



(aa 



(f){x) 



where 



2^.n 



a{zi) 



(34) 



(35) 



Note that equation (jH^ also imposes a relationship between fiix) and /j?(x) 



8 



5.2 Radion effective action 

In braneworld theories that exhibit modification of gravity at large distances, we of- 
ten find that the radion field is a ghost (see, for example, ^3 ^]). We therefore 
need to check whether this is the case in any of our models. We do this by calcu- 
lating the radion effective action to quadratic order, following JHl (see |^ for the 
supersymmetric case). 

Let us work with the left-hand coordinate patch, which has the branes positioned 
at 2; = and z = I + 6l{x), where SI = /r — /l- The metric perturbation is given by 



X 



dy I qQ'. - 
0. \y) « 



Sg^iz = Sgz 







(36) 



In order to integrate out the extra dimension, it is convenient to have both branes 
fixed. We can do this with the following coordinate transformation 



z ^ z — B{z)5l, ^ + a^{z)d^{5l) 



B{y) 



dy 



(37) 



where B{z) is some different iable function for < z < /, satisfying -B(O) = and 
B{1) = 1. While this transformation ensures that 6gf^z is still zero, the price we pay 
for fixed branes is that we now have non- vanishing Sqzz- More precisely. 



V = Xtf^ + 2- [fi + 5lB{z)] -g,, - 2a'V,VjL 



dy 

O'^iy) 



B{y) 





dy 



Sgz 



hzz = 25lB'(z) 



To quadratic order, the effective action is given by 



S^ff = -M^ I dz I d^x^ h'^^'SEab - ^ I d^x^ h'^'Sef^ 

Jo J -^^j^ Jz=z, 



(38) 
(39) 

(40) 



where 6Eab and (56/ij> are the expansions, to linear order, of the bulk and boundary 
equations of motion respectively. 



E, 



ab 



Rab — -^Rgab + ^gab 



K^2 - K^'^g^:l + 



g^ + n 



(41) 
(42) 



We find that 6Efj_u and 6Q^^l are identically zero |18^. After integrating out a total 
derivative in z, we arrive at the following effective action for the radion 



5-, 



eff 



A 



yia{zi)a'{zi) 
where AQi = Qr- Ql- 



dy 



{aa') 



d^x^g V' + 4kA' U (43) 
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5.3 "No ghost" bounds 



The radion field, 0, becomes a ghost when the coefficient of its kinetic term becomes 
negative. This would indicate a sickness in our theory. We therefore require that 



Via{zi)a'{zi)) '^^ Jq (aa'Y 



A(rT7rT37rT)--A^ / T^^O (44) 



We shall not consider what this means for de Sitter branes, because the radion is 
tachyonic, even if it is not a ghost. However, for fiat and anti-de Sitter branes it is 
worth considering this condition in more detail. 

When A = K = 0, there is no radion, and therefore nothing to worry about. When 
A = — 6A;^, K = 0, we have a(z) = e"*^'^' and the condition (1441) demands that 



1 , , 

— > — 45 

1 - 2krR 1 + 2krL 

Now consider anti-de Sitter branes with A = —Gk"^, k = —1. Given that the warp 
factor a{z) = ^ cosh [k {\z\ — c)], we can show that the condition requires that 

sinh(2fcc) sinh[2fc(;-c)] _ „;„|,foL7N 

l+2krRtaBh[k{l-c)] l+2fcrz, tanh(fcc) anLiL\^z,i\,i j ^ 

sinh[2fc(/ -c)]sinh(2/tc) ^ ^ ' 

Here we will mainly be interested in the "symmetric" case (c = 1/2), so that (j46|) 
simplifies to 

^ ^ ^2cosh(A;/) (47) 



1 + 2krR tanh{kl /2) 1 + 2krL tanh{kl /2) 

To sum up, we have shown that we can avoid ghost-like radions by placing certain 
bounds on the parameters in our model. While this is of no real use for tachyonic de 
Sitter branes, it will be important when considering fiat and anti-de Sitter branes. 



6 Gravitons 

In this section we will consider graviton perturbations in detail. We will ask two 
important questions: (i) when does the graviton become a ghost, and (ii) do we have 
a bigravity model? As with the radion, we can derive "no ghost" bounds for the 
graviton by calculating the effective action. To establish whether or not we have 
bigravity, we need to find the graviton mass spectrum. For bigravity, we expect to 
see a massless mode and an ultralight mode [SI IH El El El Q El HDl IHj • This will lead 
to modifications of gravity at large distances. However, it is important to note that 
this modified behaviour may never be observable. For example, for AdS branes there 
is a danger that the Compton wavelength of the ultralight graviton is larger than the 
horizon size [5]. Once we have established the existance of an ultralight graviton, we 
will check that it could one day lead to interesting observations. 
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6.1 Graviton mass spectrum 

Consider the metric perturbation decomposed into plane waves 

5g,,{x,z) = Y,X^;:\^)un.{z) (48) 

m 

where 

[y'-2K\')x^-;^=m\^;::;\ (49) 

From a four-dimensional perspective, Xa^i^^ corresponds to a transverse-tracefree, spin 
2 particle of mass, mjHiSllllllSlliniiZlllHliF, 50, 51j. 

We now insert ()48|) into the equations of motion (fT^ and ()24|). Since we are 
interested in plane waves, and the radion has already been accounted for, we set 
fi = T^lu = 0. The equations of motion become 

^^ + ^ + ^jM^) = 0, (50) 

in the bulk, with boundary conditions 



H T H ^ [m — 2kX 



2 9z 12M3 2a2 
at 2; = Zj. Different eigenstates are orthogonal, so that for m ^ n 
' 1 



Um{z) = (51) 



dz — - [1 + 2rL6{z) + 2rR6{z - I)] Um{z)un{z) = (52) 
I a [z) 

We can always find a zero mode that satisfies the boundary conditions 

uo{z) = Aoa^z) (53) 

where Aq is some normalisation constant. We shall now look for massive modes for 
the flat and anti-de Sitter branes of primary interest, taking each case in turn. 

6.1.1 A = K = 

Since a{z) = 1, the bulk equation is particularly simple. In addition to the zero 
mode (j53|) . we find a number of massive modes 

Um{z) = Am cos(m|2;|) + Bm sin(m|2;|) (54) 

for some constants A^, B^, that are yet to be determined. The boundary conditions 
require that 

rimAm + Bm = (55) 
[rjimcos{ml) + sin(m/)] A^ + [r/jmsin(m/) — cos(m/)] B^ = (56) 
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Note that we have used the fact that cij = 0. If we want a non-trivial solution for 
and Bm, we are only allowed quantised values of m, satisfying 

(m^rLrji — l) sin(m/) = m {r^ + tr) cos(m/). (57) 

The main KK tower will therefore be made up of states with mhoavy ^ 1/^- Are there 
any ultralight states? To answer this, consider (jFfj) when m <^l/l. We find a mode 
with mass 

_ / 6{l + rL + rji) 
"""^'^ ^ V P + mrL + rn)+QlrLrn ^''^ 

This is indeed much smaller than l/l, if we choose x = r^/l and y = tr/I so that 

1 + 3(x + y) + 6a;?/ ^ ^ 

where |e| -C 1. 

6.1.2 A = -6P; K = 

We now consider the Randall- Sundrum I model with brane-localised curvature terms. 
This time we have a{z) = e"'^'^'. By solving the bulk equation (jnOI), we again find a 
number of massive modes E21 

Um{z) = AmJ2ime''^'\/k) + BmY2ime'\'\/k) (60) 

where J„ and Yn are Bessel's functions of integer order n. The boundary conditions 
fl^ imply 

[-eiJi{me'''^/k)+mr,e'''^J2{me'''^/k)] 

+ Brn [-e^Yi{me'''^/k) + mrie'''^Y2{me'''^ / k)] = (61) 

On this occasion we have used the fact —k6i = ai/12M^. Equation ()6H) is of course 
two equations, one for i = L and one for i = R. For these to have non-trivial solutions 
for Am and B^,, we again find that m is quantised accordingly 

[Ji{m/k) + mrLJ2{m/k)] [-Yi{me^^ / k) + mr Re'^^ime^^ / k)] 

= [Yi{m/k) + mrLY2{m/k)] [-Ji{me''^ /k) + mrne^^ J2{me^^ /k)] (62) 

We now focus on the case when the brane separation is large, I >l/k. In this instance, 
the main KK tower is made up of states with mass mheavy ^ /ce~*^' jS21 IS21 EHl- Is 
there an ultralight state? Consider (jU^ when m -C ke~^\ and use the fact that 

J,(t)^i-^ + 0(t^) + (63) 



J2{t) = ^ + O {t') Y2(t) = -^2~\^ (^'l^W) (64) 
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for |t| <^ 1. To leading order, we find a mode with mass 



[kvL + \) {krn + i) e^''^ - k{l + + tr) - [kvL - \) [kvR - \) e-^^^ 

(65) 

This is much smaller than ke ^\ as long as we choose x = kri and y = krR, so that 



(x +\){y + \) e2^' - (A;/ + X + - (x - i) - i) e" 



where |e| -C 1. 

Now consider the opposite regime, when the brane separation is small, / ^ 1/k. 
The mass spectrum is well approximated by the expressions given in section 16.1.11 
when A = K = 0: the heavy modes have mass mheavy ^ 1/^, and the mass of the light 
mode is given by (|^. Note that reduces to (j3Hj) in this limit, and is therefore 
valid for all values of /. 



6.1.3 A = -6P; K = -I 

Finally, we consider AdS branes with brane-localised curvature. We shall focus on 
the "symmetric" case, where a{z) = f cosh [k (^\z\ — |)] and \ = k/ cosh [kl/ 2]. A 
detailed analysis can be found in appendix |X| There we see that the eigenstates, 
Um{z), are given in terms of hypergeo metric functions (see equation (jHH)). Again, the 
boundary conditions impose a quantisation condition on m (see equation (jlOlj) ). 

When the branes are far apart (/ > 1/k), the situation is very similar to the case 
when Tj = We have bigravity, but the Compton wavelength of the ultralight 
graviton is larger than the size of the AdS horizon. Large distance modifications of 
gravity will therefore never be observable. 

Now consider the opposite regime, when the branes are very close together (/ -C 
1/k). Again, the expressions given in section IF). 1.11 give a good estimate for the mass 
spectrum. The horizon size on either brane is given by 1/A ~ 1/k. Large distance 
modifications of gravity could therefore be observable, as long as, 

7 > "flight > k (67) 



6.2 Graviton effective action 

Now that we know the graviton mass spectrum, we are ready to calculate the cor- 
responding effective action to check for ghosts. Since we have assumed that = 0, 
we can immediately insert the metric (j48p into equation (j4(Jj) . Making use of the 
orthogonality condition we find that the effective action, to quadratic order, is 
given by 
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where 



M3 



/ dz 


Umiz) 


Jo 


a{z) 



[l + rLS{z)+rn6{z-l)] 



(69) 



As an aside, we will also examine the coupling of each mode to matter on the branes. 
To do this, we allow for T^^^-* ^ 0, giving additional terms 



Sr, 



EE 

i=L,R m 



(70) 



If we define the canonically normalised fields as 

C;irri) _ / (m) /71 \ 

we see that the mode of mass, m has the following coupling to matter on dAii, 



9. 



(0 



2 \/Cri 



(72) 



We can relate these couplings to the four- dimensional Planck mass, Mj on dAii. 
At the intermediate energy scale, mheavy 3> p ^ ^ught, gravity is due to both the 
massless and ultralight mode. We deduce that. 



96 



+ 2 



(0 

S'light 



(73) 



6.3 "No ghost" bounds 



We can now ask whether or not we have any ghost-like gravitons. If we assume that 
xt^^ is real, it is clear that this mode is not a ghost, provided 



C™^ > 



(74) 



Since we are not too concerned about the ultra-violet behaviour of this theory, we 
only need to impose this condition on the massless and ultralight modes. For the 
massless mode ()^ . this is equivalent to 



f a^iz) dz+ ria^{zi) ^ 

i=L,R 



For each case of interest, this gives 



1 



e^'^' + fcr^ - - ^ 



for A = K = 
for A = — n 



(75) 

(76) 
(77) 



k{rL + Tfj) + tanh 



kl 



kl ^ 2 
+ — sech 



^ for "symmetric" AdS branes (78) 



In addition we also require that Cught ^ 0. However, it will not be necessary to 
evaluate this explicitly in all CclSGS, clS most bigravity models will be ruled out by the 
radion bounds (jl^ and the massless mode condition (f7H|l . 
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Figure 1: Analysis for A = k = 0. 

7 Analysis 

The aim of this paper is to find a braneworld bigravity theory that is free of ghosts 
and tachyons. We saw in section El that tachyons appear on de Sitter branes in the 
form of the radion field, so we ruled these models out. Other types of brane have not 
been ruled out, and we showed in section lUTTl how one can choose parameters so that 
we get a model of bigravity. In this section we will check whether or not these choices 
are compatible with the "no ghost" bounds of sections 1^31 and We will also keep 
in mind the validity bound 0, and the hope that bigravity will lead to modifications 
of gravity that will one day be observable. 

7.1 A = K = 

Let us begin with the case A = k, = 0. For bigravity, we need equation (|59|) to hold 
for |e| ^ 1 (recall that x = r^/l and y = r^/l). We have no radion to worry about, 
although we do need to be careful with the gravitons. The massless graviton is not a 
ghost, as long as 

l + x + y^Q (79) 

Now consider the plot oi x = ri/l against y = r^/l shown in figure ^ The solid black 
line corresponds to a line of constant |e|-Cl. Asx— >-oo,?/— l/6e^, and vice- versa. 
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The dotted line corresponds to 1 + x + y = 0, and represents the boundary between 
the massless graviton being well behaved, and the appearance of a ghost. Since the 
solid line never crosses the dotted line, it is clear that we need x, y l/6e^. In other 
words, we need rj 3> /. 

We are now ready to check the behaviour of the ultralight mode. For simplicity, 
let = = r /, so that 

flight ~ (80) 

If we use this along with equations ()54p and ()55p . we find that 

Qight ~ AlM^r ^ (81) 

The ultralight mode is not ghost, so we conclude that we have indeed found a bigravity 
model that is ghost-free. At intermediate energy scales, note that the Planck mass 
on either brane is given by 

Mf = mil ^ (82) 

as we might have expected. 



7.2 A= -6A;2; / <Cl//c; «: = 0, -1 

We now turn our attention to the case where A = — k = 0, —1, with I <^ 1/k. In 
this limit, the mass spectrum is the same as for A = k = 0, and the graviton bounds 
(f77|) and (fTHj) . are reduced to (fT^ . The analysis will therefore be very similar to 
what we have just done for A = k = 0, as we would expect by continuity. There are, 
however, a couple of other things to consider now. We have a radion to worry about 
(see equations and (jlZj)), and for k = —1, we need to ensure that any infra-red 
modifications of gravity will not be hidden behind the AdS horizon. 

For simplicity, let us again assume that = rji = r ^ I. In both cases we find 
that there is an upper bound on the value of r 



r< , (83) 




For K = 0, this is due to the radion becoming a ghost jSH!- For k = —1, however, 
the radion is not a problem, and the bound is due to the condition, mught > k. This 
ensures that long distance modifiactions of gravity may one day be observable. 

The upper bounds on r are not too troublesome, as long as we take k to be very 
small. 



7.3 A=-6P; l>l/k; k = 

Finally, we consider the (DGP) extension of the RS model, with the branes far apart. 
For bigravity, we need equation to hold for |e| ^ 1 (recall that now we have 
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, No 
ghosts 



Both ghosts 



l^adion ghost 
only\ 



1/2 




Graviton ghost 
only 



No ghosts 



Figure 2: Analysis for A = —Gk"^, k = 0, I >l/k. 



X = kri and y = kr^. From equations (jl^ and (f77j) we have the following bounds 

-2kl 



-1 e 
> — 



y 2 + 2 



to avoid a radion ghost 



(84) 



— - )e^*^'^0 to avoid a (massless) graviton ghost (85) 



Now consider the plot oi x = kr^ against y = kr^f shown in figure |21 Again, the solid 
black line corresponds to a line of constant |e| <^ 1. The asymptotic behaviour is 



X 



1 + coth(2A;/) 

cosech(2A;/) 
2? 



as X ^ oo 



as 1/ — >• oo 



The dotted lines correspond to boundaries, across which the "ghost-like" status of 
the theory changes. The diagonal hne is given bya; + i + {y ~ \) e"^'^' = 0, whereas 
the horizontal and vertical lines are given by ?/ = 1/2 and x = —1/2 respectively. It is 
easy to check that the solid line never enters any of the "no ghost" regions indicated. 
This means that for this model to exhibit bigravity, either the radion, the graviton, 
or both, will become a ghost. The model is therefore rejected. 
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8 Conclusions 



In this paper we have studied a class of braneworld models, with localised curvature 
on the branes. Some of these models give rise to bigravity, leading to large distance 
modifications of gravity for a four- dimensional observer. After checking for ghosts 
and tachyons, we rejected most, but not all of the models we had considered. 

One of these well behaved models consists of a flat bulk (A = 0) sandwiched in 
between two flat branes {k = 0) of zero tension. The branes are close together, and 
the curvature terms on the brane are large and positive (r^ ^ /). The mass spectrum 
is precisely that of bigravity. We have a massless graviton, an ultralight graviton, and 
a tower of heavy KK modes. The model is completely free of ghosts and tachyons. 

To add some substance, let us put some numbers in. If we take I ~ 10^ (eV)~^, 
and Tj ~ 10^^ (eV)""^, we find that 

miight ~ 10-3^ eV, mheavy > 10-' eV. (86) 

These masses lie outside of the range for which gravity is well tested (H)), so we 
have no contradiction with experiment. Since the four-dimensional Planck mass, 
rupi ~ 10^^ GeV, we conclude from equation (jHH) that M ~ 10"^ - 10"^ eV. While 
the fundamental Planck scale is very low, it does not violate the validity bound 0, 
and exceeds the scale probed by Cavendish experiments (10~^ eV). In principle there 
are cosmological and astrophysical bounds that one should consider, such as the effect 
on star cooling due to graviton emmision into the bulk. This will be left for future 
research. 

Note that mugut ~ 10~^^ eV is just about heavy enough for us to expect interesting 
observations today. This could be important in trying to explain cosmic acceleration. 
If we wanted to reduce this mass even further, we would have to increase r. Even a 
small increase would push the Planck mass below the Cavendish scale. Curiously, we 
seem to be in just the right place to start seeing interesting new physics, either in the 
infra-red, or the ultra- violet^. 

We can go beyond this model by switching on a small AdS curvature in the bulk, 
and still get ghost-free bigravity. Although a tachyon appears for de Sitter branes, 
this is not the case for flat and anti-de Sitter branes. When the bulk AdS length is 
of the same order as the brane separation, we find that ghost-free bigravity becomes 
impossible. 

By checking for ghosts and tachyons in our models, we have carried out the first 
important tests of viability. However, we should be aware that there are other possible 
problems. As with all models of massive gravity, we should be concerned with the 
famous vDVZ discontinuity |23I2H|- Our linearised equations of motion suggest 
that this may well be an issue. This is especially true for the case A = k = 0, as there 
is no "brane-bending" effect that could cancel off any unwanted degrees of freedom 
(see, for example, ESI)- Even our AdS brane model will suffer this problem. 
This is a surprise, as it is often said that there is no vDVZ discontinuity in AdS 
space jSniEZI- However, this result relies on the graviton mass going to zero faster 

would like to thank John March-Russell for this observation 
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than the inverse horizon size. In order to ensure that we had observable bigravity, we 
demanded the opposite of this. 

Crucially, there may be a way around this problem. Vainshtein et al |^ IH^ . 
have argued that the perturbative expansion in Newton's constant in [23 is 
inconsistent, as the graviton mass goes to zero. Specifically, the standard linearised 
analysis near a heavy source is only valid at distances 



where tm is the Schwarzschild radius of the source, and m is the graviton mass. To 
see if a mass discontinuity really is present, we need to calculate the Schwarzschild 
solution on a brane, and compare it to the standard four-dimensional massless result. 
This is a highly non-trivial exercise that is beyond the scope of this paper. 

For Pauli-Fierz theory , the breakdown of the linearised analysis has been linked 
to strong coupling phenomena at the following scale pTj 



where Ipi is the Planck length^. Large terms in the full propagator tend to signal this 
problem. Unfortunately, we might expect our theory to suffer the same fate. We can 
see this by looking at h^z in a fixed wall gauge (see equation (jHOl))- By the mean 
value theorem, there exists zq G [0,/] such that B'{zo) = l/l. For small I, it is clear 
that hzz can be very large. 

A similar strong coupling scale also exists in the DGP model jHUl |M] . However, it 
has recently been argued that this scale could be unphysical, and is just a result of the 
naive perturbative expansion jHZI- Clearly, both the mass discontiniuty and strong 
coupling problems are highly contentious issues at the moment. For this reason we 
have focused on the possible existence of ghosts in our models. There is no contention 
there: ghosts are undesirable, but can be avoided. 

There is still much to do. It would be very interesting to study the phenomenology 
of these models in more detail, particularly in the context of cosmic acceleration. 
Coincidentally, we have been forced to choose a brane separation / ~ 10'^ (eV)~^, 
which agrees with the brane separation chosen in In the mass spectrum is 
motivated by a discretized Randall-Sundrum model, with / chosen so that there is a 
small effective cosmo logical constant. Would the presence of an ultralight mode in 
the spectrum affect these results? We would probably expect the answer to be "no" , 
because the dominant contribution to the vacuum energy would still come from the 
heavier modes. Nevertheless, it is worthy of further investigation. 

Our models could certainly be generalised in a number of ways. We could consider 
more branes, abandon Z2 symmetry, or even introduce higher derivative terms in the 
bulk and on the brane. Branes embedded in solutions to Gauss-Bonnet gravity have 

■^For certain non-linear extensions of Pauli-Fierz, i?strong ~ (V/'^^) " i2fi|- This higher scale 
also appears in the DGP model |6U|. 




(87) 




8. 
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been the subject of much research recently (see, for example j64[ I65j). motivated by 
the link to string theory. Indeed, it would also be nice if our models, or at least some 
generalisation, could be derived from a more fundamental theory [221 12n] • However, 
the high degree of fine tuning could prove an obstacle in this respect. 

Let us end by summarising our main result: we have discovered braneworld models 
that exhibit bigravity, without introducing ghosts. Recall that bigravity naturally 
gives rise to new gravitational physics in the infra-red. The new physics occurs when 
the massive graviton "switches off", at distances beyond its Compton wavelength. 
Our models are an improvement on the ghost-free model given in ^j, because they 
lead to potentially observable modifications of gravity. In P, all modifications are 
hidden behind the AdS horizon. 
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A Mass spectrum for AdS branes 

For "symmetric" AdS branes, we have 



We want to find the massive eigenstates, Um{z), that satisfy equation (jSOl)- 




A = A;/ cosh [kl/2] 



(89) 




(90) 
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Let Um{z) = Vm{y) /ci^iz), wheie y = tanh^ (A; [\z\ — |)). We find that Vm satisfies 
the following hypergeometric equation 



1 - AX^ 

1/(1 - yWLiy) + - '^y>Liy) + -^^^Uy) = 

The solution is given in terms of hj^ergeometric functions jO] 
where Am, Bm are arbitrary constants and 



The eigenstates are therefore given by 



Um{z) = sech^ 
+Bm tanh 
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where we have absorbed some constants into Am and Bm- After some tedious algebra, 
we see that the boundary conditions (j51|) give the following 



AmO^m ~l~ Bmf^m 
AmO^m Bm(3m 



(95) 
(96) 



where 



Or 



Aktambm {t^ -l) F (^m + 1, 6m + 1, t' 

+ {Akt + fm^) F (^m, bm, ^; 



m 



Ar 



1 3 



a 



m 



—j m tF l^ttm + -,bm + -,-;t 

n^f , 1 ^ 

— ]m F [ am,bm, -]t' 



(97) 
(98) 
(99) 



m 



^kt" [e - 1) (^m + {b^ + F (^a^ + ^, 6„ + ^, ^; 
+ (A; (5^2 - 1) + fmH) F (^m + ^, + ^ | 



(100) 



Here t = tanh(A;//2) and r = (r^, + rR)/2. By demanding that there exists a non- 
trivial solution for Am and i?„i, we find a quantisation condition for m 



(101) 
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